Theoretical and computational study of high pressure structures in barium 



Stewart K. Reed and Graeme J. Ackland 
Department of Physics and Astronomy, The University of Edinburgh 
Edinburgh, EH9 3JZ, Scotland, United Kingdom 

Recent high pressure work has suggested that elemental barium forms a high pressure self-hosting 
structure (Ba IV) involving two 'types' of barium atom. Uniquely among reported elemental struc- 
tures it cannot be described by a single crystalline lattice, instead involving two interpenetrating 
incommensurate lattices. In this letter we report pseudopotential calculations demonstrating the 
stability and the potentially disordered nature of the 'guest' structure. Using band structures and 
nearly-free electron theory we relate the appearance of Ba IV to an instability in the close-packed 
structure, demonstrate that it has a zero energy vibrational mode, and speculate about the struc- 
ture's stability in other divalent elements. 
PACS numbers: 71.15.Hx, 61.50.K, 71.20.Dg 



The high pressure structures of the alkaline earth met- 
als show an unusual trend in going from apparently close 
packed structures at low pressure to more complex open 
structures at higher pressure. This counterintuitive be- 
havior arises from transfer of electrons from free-electron 
like s-bands to more directional d-bands. 

Barium has a number of still more unusual features. 
It adopts the body-centered-cubic (bee) structure at 
zero pressure (conventionally denoted phase I). Under 
5.5 GPa of pressure it transforms to phase II with a 
hexagonal-close-packed structure (hep). At 12.6 GPa it 
transforms into the complex phase IV and then back to 
hep (phase V) at 45 GPa. This last phase is observed 
experimentally up to at least 90 GPa. Until recently all 
attempts to solve the phase IV structure have failed. 

Nelmes et al pj have now shown that Ba IV has a 
complicated structure consisting of two interpenetrating 
but incommensurate structures: a tetragonal 'host' with 
'guest' chains in channels along the c-axis of the host. 
Although they share the same a and b parameters the 
guest and host lattices are incommensurate in the c-axis. 
Viewed along c, the eight atoms in the host unit cell are 
arranged in an octagon but with atoms at alternating 
heights. The guest chains pass through the middle of 
these octagons, and a combination of diffuse scattering 
and Bragg peaks suggests regions of both order and disor- 
der between the chains. Figure [l] illustrates the unit cells 
and describes how Ba IV splits up in distinct 'subphases' 
depending on the registry of adjacent chains. 

An unusual feature of the hep phase is the strong pres- 
sure dependence shown by the c/a ratio in phase II: it 
falls to ^1.50 at the II-IV transition, so far below the 
ideal value (1.633) that the structure can hardly be re- 
garded as close packed. In phase V however, which is 
also hep, the reported c/a ratio is almost pressure inde- 
pendent at ~1.58 0. 

In this paper we present total energy calculations on 
the observed and other possible structures. These en- 
able us to probe experimentally unobservable quantities 
such as the energy barrier to moving the chains and 



whether there is an isostructural II- V phase transition 
'hidden' in the range of pressures where Ba IV is more 
stable. Dealing with the incommensurate double-lattice 
provides a unique challenge for simulations based on pe- 
riodic boundary conditions. 




FIG. 1. Schematic representation of the Ba IV structure. 
The host atoms (dark symbols) with chains of guest atoms 
(light symbols) are projected on to the ab-plane. The host 
atoms are labeled with their z coordinates. The arrows la- 
belled a„ and b s are projections of the supercell vectors used 
in the calculations. At 12.1 GPa the host cell has dimensions 
o-host=bhost=8-42A and Chost=4-74A. Full description of the 
observed guest atoms involves three unit cells. One is par- 
tially disordered and tetragonal with the same a and b as the 
host structure but with c gueB t t =3.J^lA. On ordering, the guest 
structure can undergo a monoclinic distortion to one of four 
similar cells, related by 90° rotations about the tetragonal 
axis, with a guestm =8.46A, b guestm =a host , c gues t m =3.4SA 
and f3 guestm =96.15° . This corresponds to displacement of 
atoms in adjacent chains by O.J^lA along c. 

The calculations are done using the ab initio total- 
energy density- functional-theory |p -pfl plane- wave ultra- 
soft pseudopotential method |||l0 which is well docu- 
mented elsewhere. The generalized gradient approxima- 
tion [Q was used for exchange and correlation. The plane 
wave basis is cut off at 350 eV which is sufficient to con- 
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verge total energies to better than 0.1 meV per atom. 
Sampling of the reciprocal space uses symmetry-reduced 
k-point sets H converged to 0.1 meV per atom. Eight 
valence electrons per atom are treated explicitly. We cal- 
culate the total energies for the bcc, fee, hep and Ba IV 
structures. In all cases the enthalpy is minimized with 
respect to ionic and unit cell degrees of freedom under 
hydrostatic pressure using Hellman-Feynman forces and 
analytic stresses with Pulay corrections ||. 

Plane wave calculations are done with periodic 
boundary conditions, but it is impossible to model 
an incommensurate structure exactly using a super- 
cell. Hence Ba IV was approximated by a simi- 
lar commensurate structure: pseudo-Ba IV comprises 
a monoclinic supercell (s) containing eight host and 
three guest atoms with a s = \{b host + a host ) + ^c hosU 
b s = \{bhost - ahost) + \chost and c s = 2c host . Before re- 
laxation three guest atoms were evenly distributed along 
the c-axis from the origin. This arrangement gives a Chost 
to Cg Ues t ratio of 1.5 compared to experimental value of 
about 1.39. Furthermore, since a s and b s both contain a 
component in the Chost direction the neighboring columns 
of guest atoms are displaced in the Chost direction with 
resp ect to each other more than is seen experimentally 
jllj . The total energy for this cell will represent an upper 
bound on the energy of incommensurate Ba IV. 
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FIG. 2. Graph of energy against volume for various 
phases. At larger volumes bcc has the lowest energy, while 
hep is lower below ~ 48 ± 2 A 3 . The inset shows the differ- 
ences in enthalpy between the hep, BalV and fee with refer- 
ence to that of bcc. The calculated Ba IV data is shifted by 
-0.065eV (see text) to account for the incommensurability. 



The results of the calculations are shown in hgure |^. 
The energy differences between the bcc, fee and hep 
structures are small but they are nonetheless greater than 
the accuracy of the calculations. In agreement with ex- 
periment, bcc is the favored structure at ambient pres- 
sure with a phase transition to hep at around 5 GPa. 



The pseudo-Ba IV phase energy is close to the hep curve, 
moving away at low and high pressure: evidence of the 
reentrant nature of the hep phase. 

The overestimate of the energy in calculating pseudo- 
Ba IV rather than incommensurate Ba IV can be approx- 
imated by a simple elastic model: assume that to form 
pseudo-Ba IV the Ba IV guest lattice has been indepen- 
dently compressed and the host lattice expanded until 
the experimentally measured ratio Ch ost / c guest =1.39 0] 
reaches the commensurate value of the supercell: 1.5. 

By measuring the stress 033 of the commensurate cell 
under uniaxial compression £33 and assuming that guest 
and host lattices separately have the same elastic con- 
stant C33, linear elasticity yields a strain energy of: 
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The composite elastic modulus C33 can be measured 
by a finite strain Jl2] calculation, taking a pressure in- 
dependent C33=150GPa gives the energy shift AE = 
0.065eV/atom which is included in figure 2. The Ba IV 
enthalpy lies below hep between 9.5 GPa and 23 GPa, in 
reasonable agreement with experiment. The exact tran- 
sition pressures are sensitive to this correction. 




FIG. 3. In an incommensurate structure, guest atoms oc- 
cur at all possible heights in the host unit cell. The guest-host 
interactions could induce a modulation with period c gues t in 
the host which would yield an indirect chain-chain coupling 
(upper figure, with the host displacements much exaggerated). 
Since all local registries are observed somewhere along the 
chain, the energy change due to a vibrational mode displacing 
the guest chains along their length with respect to the host is 
exactly zero (lower figure) . The host modulation would move 
with the guest atoms holding adjacent chains in registry with 
one another. 



The offset of the origin of the two lattices must be 
defined for commensurate structures: we used a range 
of origin offsets and found that the energies varied be- 
tween them by less than O.OOleV/atom and the relaxed 
volumes by less than 0.02% |ll| . This shows that no posi- 
tion of the guests is more favorable than any other which 
means that the guests can be at any height with respect 
to the hosts without there being an 'energy penalty' to 
pay. Thus the c gues t is determined by the packing along 
the chain rather than host-guest interactions. In the fully 
incommensurate case this would lead to a zero frequency 
phonon mode, or phason (figure 
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The calculated c/a ratio was 0.58 as compared with 
experiment Ch os t/o,host value of 0.56. This can be ex- 
plained by the elastic model: 0.2 more guest atoms have 
been forced into the cell than would be ideal, applying a 
stress on the host. The host cell therefore responds by 
expanding in the c-direction. 

Moreover, hep calculations show a plateau in the Fermi 
energy between 18 GPa and 43 GPa resulting from a peak 
in the density of states. Typically, this would lead to a 
Jahn- Teller type instability of hep. 

Recent studies [^],|l4j of the two hep phases have shown 
a dramatic decrease in the c/a ratio in phase II with in- 
creasing pressure. In phase V it is noted that the ratio is 
almost constant, slightly less than ideal. Experimentally 
it is impossible to get data between 12 GPa and 45 GPa 
because hep is unstable in this region; it is nevertheless 
possible to do calculations. 
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FIG. 4. Pressure dependence of the Fermi energy (top) 
and c/a ratio (bottom) for hep Ba, from constant (hydro- 
static) pressure calculations. The open triangles are experi- 
mental data taken by Takemura jj/. 



As figure |j illustrates our c/a results initially agree 
very well with experiments and remain in good agree- 
ment up to the phase transition. In the experimentally 
inaccessible region we find a steady increase in c/a, and 
decrease in volume until once again good agreement is ob- 
tained when the hep phase is again observed at 45GPa. 
Thus we have demonstrated that phase II and phase V 
are actually the same phase, with the variation in c/a due 
to a continuous transfer of electrons from s to d. Zeng et 
al calculated the c/a ratio using the linear muffin-tin 
orbital method for electrons in a variety of different con- 
figurations. They found that, at 5.7 GPa near the onset 
of phase II, a better fit to the experiments was obtained 
if they included d-orbitals as well as s- and p-orbitals. 
They attribute the variable c/a ratio to strong s — > d 
transfer across phase II, but report a much smaller effect 
in phase V. 

To interpret these results, we consider packing of hard 
spheres and the nearly free electron model. (l6). As the 
host and guest cells are incommensurate the guest atoms 
must be able to take take any (every) height with respect 



to the host cell. In the packing model this implies that 
the guest atoms must be able to pass through the octago- 
nal host rings. With hard sphere atoms this would mean 
either the host atoms cannot be in mutual contact or the 
guest atoms must smaller than the hosts. Even then, the 
packing fraction is 0.56 whereas, at the II-IV transition 
the hep phase II has a packing fraction of 0.68. This sug- 
gests that efficient packing is insufficient to understand 
the structure. 

In the electron gas/pseudopotential perturbation the- 
ory model of Heine and Hafner jl6| , screening of the nu- 
clear charges gives rise to an oscillatory effective pair po- 
tential between atoms. The favored crystal structure is 
determined by this pair potential, and for Ba the high 
pressure close packed structures have atoms at unfavor- 
able separations, near maxima in the pair potential p5[ . 
Thus the 'ideal' hep structure is unstable with respect to 
splitting the shell of 12 nearest neighbors into 6 and 6 by 
altering the c/a ratio. 
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FIG. 5. Observed bond lengths as a function of pressure 
in hep, bec and Ba IV structures. The Ba IV interhost dis- 
tances are similar to the neighbor distances in hep, which 
themselves are a continuation of the bec distances. The Ba 
IV guest separation is consistently lower. The shortest pos- 
sible guest-host distance is half the distance across the host 
octagon, which is slightly more than the smallest guest-guest 
spacing. Hence the guest atom is more compressed than the 
host. 

Since the potential depends only on interatomic spac- 
ings, competing structures will have neighbor separations 
similar to the optimally distorted hep. In our calculations 
( figure H) we find that bec and hep have approximately 
the same nearest-neighbor distance at room temperature. 
As the pressure increases the hep nearest neighbor dis- 
tance splits into two due to the non-ideal c/a ratio. At 
the II-IV transition the intra-guest distance is 3.4lA, and 
the shortest intra-host distance is 3.43A. These distances 
are close to the nearest neighbor distances in hep (3.46A). 
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As a result of the incommensurability, a range of host- 
guest separations nearest-neighbor distances should oc- 
cur. The similarity of bondlengths (figure ||) to those of 
hep shows that Ba IV structure represents another way 
of arranging the atoms such that they lie in a minima of 
the screened electrostatic potential. 

An alternate view of the nearly free electron picture 
jL7| is that the energy is lowered by perturbation of states 
near the Fermi surface: this favors structures which have 
Brillouin zone faces close to the Fermi vector. Incom- 
mensurate Ba IV has two Brillouin zones with the (220), 
(211) host and (lll) gU est facets of the first zones close to 
the divalent free electron Fermi vector |L8| . 

In sum, the total energy pseudopotential method gives 
an excellent description of the experimentally observed 
high pressure phase diagram for barium, including the 
complex guest-host structure Ba IV. Although one can 
not represent the full incommensurate structure in a su- 
percell calculation, the energy of an equivalent commen- 
surate pseudo-Ba IV structure is close to that of hep: 
treating the mismatch using linear elasticity shows that 
incommensurate Ba IV has a region of stability. 

Calculations involving moving the guest chains within 
their channels, show the interchain interactions of the 
guest atoms to be very weak, probably mediated by elas- 
tic strain in the host. Based on this we propose that 
the tetragonal phase (lj is metastable. Weak interactions 
mediated by strain in the host lead to ordering at low T. 
The least favored arrangement in terms of strain would 
be to have all guest atoms simultaneously coplanar, and 
so each chain is displaced relative to its neighbors leading 
to the monoclinic phase |p^| . 

Ba IV an unusual example of a stable elemental struc- 
ture with two inequivalent atomic sites. Pressurization 
of barium proceeds by transfer of electrons from low en- 
ergy, extended s-states to higher energy, more localized 
d-states. The hep structure represents one method of 
achieving intermediate sd bonding: hybridized orbitals 
with all atoms equivalent. Ba IV offers another mech- 
anism, with a division between guest atoms containing 
more d-like character and host atoms having more s-like 
behavior: a projection of the Kohn-Sham wavefunctions 
onto atom centered orbitals shows s — d transfer in all 
structures and a 20% greater s-character on the host 
atoms compared to the guest atoms in the region of Ba 
IV stability. The p and d-character discrepancy was less 
pronounced, suggestive of a transfer from s to more free- 
electron like behavior. The larger number of electrons 
'localized' on the host atoms is consistent with their ap- 
parently smaller size. Thus Ba IV can be regarded as an 
intermetallic compound, in which both components are 
actually the same element. 

The free electron theory suggests that Ba IV will be a 
competitor with hep as a high pressure phase not just in 
Ba, but in other materials with similar valence electron 
density and s — > d electron transfer. Moreover the re- 
duced size of the d-like atom will lead to faster diffusion 
and lower melting points at high pressures. 
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